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<^ ■ Abstract 
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^ ! A conjectured duality between supergravity and N = oo gauge theories 

gives predictions for the glueball masses as eigenvalues for a supergravity wave 
' equations in a black hole geometry, and describes a physics, most close toa high- 

, temeperature expansion of a lattice QCD. We present an analytical solution 

Q>^ I for eigenvalues and eigenfunctions, with eigenvalues given by zeroes of a cer- 

tain well-computable function r{p), the zeroes of which signifies that two solu- 
tions with desired behaviour at two singular points become linearly dependent. 
Qh! Our computation shows corrections to the WKB formula = 6n(?7- -|- 1) for 

eigenvalues corresponding to glueball masses in 3 dimensional QCD, and gives 
the first states with masses = 11.58766; 34.52698; 68.974962 ; 114.91044; 
172.33171; 241.236607; 321.626549, .... In QCD^, our computation gives 
X ; squares of masses 37.169908; 81.354363; 138.473573; 208.859215; 292.583628; 

^ '• 389.671368 500.132850, 623.97315 ... for 0++. In both cases, we have a power- 

ful method which allows to compute eigenvalues with an arbitrary precision, if 
so needed, which may provide quantative tests for the duality conjecture. Our 
results matches well with the numerical computation of ||5| withing precision 
reported there in both QCD3 and QCD4 cases. As an additional curiosity, we 
report that for eigenvalues of about 7000, the power series, although conver- 
gent, has coefficients of orders lO^'*; and also the final answer gets small, of 
order 10"*^ in QCD4. Tricks were used to get reliably the function r{p) above 
7000. In principle we can go to infinitely high eigenavalues at an expence of 
computer sufferings; although eventually such computation will slow down to 
make it inpractical, also as corrections may be expected for higher states, since 
fine cancellations of very big terms is what produces higher eigenvalues. 
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1 Glueball masses in QCDs 



The conjectures of Maldacena and Witten 0|, [|T] gives, in particular, a prediction 
for glueball masses in iV = oo QCD in 3 dimensions from supergravity, namely, from 
certain solutions of classical equations of motion in a black hole metric for a massless 
scalar dilaton field $. The dilaton shows up here as it couples to the 0'^'^ glueball 
operator; other glueballs, which couples to other fields, say a two-form in supergravity 
theory, are described, for example, in . The classical equations of motion for dilaton 
are: 



in a black hole metric 
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Here Xi,X2-,x^ is where our QCD^ lives, and it roughly should be imagined as a 
4-ball Xi^X2tX^, p with xi^X2,x^ on boundary; r is a circle, on which we pose 
antiperiodic boundary conditions for fermions and periodic for bosons, thus breaking 
supersymmetry; dVt^ is a sphere with a standart metric. We need to take solutions 
of the form 

with the appropriate physical conditions for $ at p = oo, (normalizability) , and at 
p = 1 (single- valuedness for complex p) By introducing a variable X = p^, (and 
calling it just x below) we get certain well-posed Sturm-Liouville problem: 



dx"^ \x x — 1 x + lj dx x{x 
together with the boundary data: 

/ decays as x ^ oo, and: 

/ is regular at x = 1. Here Ap = fc^ is the mass of glueball states, Ap = —k'^. 

The above equation has regular singularities, and its solution can be written as a 
series, with radius of convergence determined by distance between two singularities. 
For the above equation, the radiuses of convergencies of expansions at 1 and infinity 
overlap, therefore, a reliable computation is possible. 

For an arbitrary p, the fundamental system of solutions of the equation ( |l|) is the 
following: 

For X with 1 < x < oo, it is given by a linear combination of two convergent series 
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oo 



a\n\x 



n=l 

oo 



(2) 



n\x 



n=l 



where a[n], b[n] are given by a recursion formula 

1 9 

a[< 0] = 0;a[0] = l;a[n + 1] = ^ ia[n - l]((n + 1) ) + pa[n]); 

{n + 2) - 1 

b[< 0] = 0; b[0] = 0; 

+ 1] = ^ (— (2na[ra- 1] - {2n - 2)a[n - 3]) + b[n - l]{n - if + pb[n]] 

(n) - 1 V 2 / 

For X with < x < 2, the fundamental system of solutions is given by a linear 
combination of two convergent series 

oo 

y['\x) = l + Y.a[n]ix-ir 

^ oo ^ (3) 

n=l 

where a[n], 6[n] are given by a recursion formula 

a[< 0] = 0; a[0] = 1; a[n + 1] = ^ ((3n(n + 11) - p)a[n] + a[n - \\[n^ - 1)) 

2{n + 1) ^ ^ 

b\< 0] = 0;6[0] = 0; 

b\n + 1] = ^ 2 (^('^ + + 1] + 3(1 + 2n)a[n] + 2na[n - 1] + 

2{n + 1) 

+ (3ra(n + 1) - v)b\n\ + 6[r2 - \\{n^ - 1)) 

The solution of the boundary problem ([1|) can be described as follows: For certain 
p there exist a solution which is proportional to y^^^{x) for x > 1 and to y^\^{x) as 
|a; — 1| < 1. The condition for such a solution to exist is that the Wronskian of two 
solutions y^^\x) and y^Y^{x\ which depend from p as a parameter, 



Wronskain{j)^ x) 



/ y\'\x) yr{x) \ 

must be zero. For x such that 1 < x < 2 both the series defining y^^\x) and the 
series defining y^i\x) are convergent, therefore, the Wronskian can be computed. It 
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is easy to see that the Wronskian depend from X as 



Wronskain{p, x) 



r{p) 



x{x — l)(a; + 1) ' 



and therefore, the function r{p) 



r{p) = x{x — l)(x + l)Wronskain{p, x) 



can be computed at any point x with 1 < x < 2, which allows to determine the 
function r{p), with any desired accuracy, as the series (H), (H) are convergent, and 
also the solutions yi are analytic functions of x. 

The spectrum of ( ??) correspond to zeroes of the function r{p). By standart 
oscillation theorems, there are no zeroes of r{p) for p > 0, and there is an infi- 
nite discrete set of zeroes for p < 0. We do not know an analytic formula for the 
roots of r{p), however, we can find their numerical values, for example by plot- 
ting function r{p) acurately and looking where the zeroes are, as one would do 
for a transcendental equation of the sort tan(x) = x; we can also find roots nu- 
merically with an arbitrary precision, for example by invoking Newton method of 
finding a root. Such computations shows that the first several roots are located at 
-11.58766; -34.52698; -68.974962; -114.91044; -172.33171; -241.236607; -321.626549,.... 
We also believe that numbers Ap = —Qn{n + 1) , n = 1, 2..., which WKB method give, 
are not the roots, although they are close to those roots listed. 

1.1 Orthogonality 

For the discrete set of values p„ , n = 1, 2 . . . such that r{p) = we have normalizable 
at 1 < X < cx) wave functions . Those functions are orthogonal, which here 



1.2 Checking that -12 is not a root 

Since a WKB computation, which may or maynot have corrections in this case, shows 
the first eigenstate at Ap = —12, we were interested to check does our function r{p) 
has a root exactly at 4p = —12. It is easy to check whether or not it does, as our 
series are convergent (and pretty fast, faster then the worst of of X](^l)"'^(3; — 1)" and 
J2{—^)"'n{l/x)'^ do), for the p we are interested in). Therefore, we can compute the 
function r{p), using finite number of terms in the series, and estimate the error. Thus, 



says 




(4) 
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using our method, we can answer the question whether or not the WKB formula gets 
corrections. We did such computation, and our result is that 

-0.022482 < r(p)|4p^_i2 < -0.022481 

(with the sign of r{p) here consistent with our expectation to get the first root close 
by with 4p > —12, as r{p) is positive and of order 1 for small negative p). Thus we 
believe there is no root at exactly 4p = —12, and there is a correction to the WKB 
formula, with the exact root at about —11.59. 



2 Glueball spectrum in QCD4 

The conjectures of Maldacena and Witten |^, |l[| gives, in particular, a prediction 
for glueball masses in = oo QCD in 4 dimensions from supergravity, namely, from 
certain solutions of classical equations of motion in a black hole metric for a massless 
scalar dilaton field $. The dilaton shows up here as it couples to the 0~^^ glueball 
operator; other glueballs, which couples to other fields, say a two-form in supergravity 
theory, are described, for example, in . The classical equations of motion for dilaton 
are: 

where the metric is a black hole 

Is 95^ N/ An \ Py t^i 

Here where our QCD4 lives, and, oversimplifing a bit, it can be imag- 

ined as a 5-ball xi, X2, X3, X4, p with xi,X2,X3,X4 on S"^ boundary; r is a circle and 
dQ^ is a sphere S'^ with a standart metric. We need to take solutions of the form 

with the appropriate physical conditions for $ at p = 00, (normalizability) , and at 
p = 1 (single- valuedness for complex p) By introducing a variable X = p^, (and 
calling it just x below) we get certain well-posed Sturm-Liouville problem: 

d'^F / ^ dF 

{x^ -x)^ + (lOx^-A)^-px^F = 0, (5) 
dx'^ ^ ^ dx 

(which is the equation of motion of a massles scalar field, dilaton, which couples to 
the relevant glueball operator) 
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together with the boundary data: 
F{x) decays as x — oo, and: 

F{x) is regular at x = 1 Similarly to the equation ( |l|), the singularities at 1 and 
infinity for this equation are regular, and the fundamental system of solutions of the 
equation ( ^ can be written as power series with the appropriate radius of conver- 
gence, namely, for 1 < x < oo the fundamental system of solutions is is 



y[°°\x) = l/(x9) + ^ a[-i--^-" 



yt\x) = i^Log{x)y[^\x) + Y.b[ 



n\x 



n=l 

oo 



(6) 



n\x 



n=l 



where a[n] are given by a recursion formula 

a[< 0] = 0;a[0] = 1; 

a[2n] = ^^^^^l^^g^ (P«[2^ - 2] + (2n)(2n + 3)a[2n - 6]) ; n = 1, 2, . . . 
The series are convergent for 1 < x < oo. 

For X with |x — 1| < ^ — , (with r= — being the distance between 1 and the 

' ' ~ V V2 

next nearest sixth roots of one, e~6~, on a complex plane of x, as can be discovered 
by looking at the coefficient in front of the second derivative in our equation) the 
fundamental system of solutions is given by a linear combination of two convergent 



series 



yi'\x) = l + Y.a[n]{x-ir 

^ 00 ^ (7) 

yi^\x) = Log{x - l)y[^\x) + ^ b[n]{x - l)"", 

n=l 

where a[n] are given by a recursion formula 

a[< 0] = 0;a[0] = 1; 

+ ^] = -7^r^ri\^ (21(n - 1) + 60) - p)a[n] + 

{{n - 1) (35(n - 2) + 150) - 3p)a[n - 1] + ( (n - 2) {35{n - 3) + 200) - 3p)a[n - 2] + 
+ ( (n - 3) (21(n - 4) + 150) - p)a[n - 3] + {{n - 4)(7(n - 5) + 60))a[n - 4] + 
+ {{n - 5){n - 6 + 10))a[n - 5]) 
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Our boundary conditions require that the solution must be proportional y[°°\x) for 
X > 1 and to as |a; — 1| < 1, and therefore, the Wronskian of two solutions 

yi^^(x-) and yi \x), 



Wronskain{p, x) = 

which depend from p as a parameter, must be zero. For x such that 1 < x < 

the series defining y[°°\x) and the scries defining y^\x) are both convergent, (for 
any p), and therefore the Wronskian can be effectively computed for such x using our 
series. The Wronskian depend from x as follows: 

r(p) 

Wronskain{p,x) — 



x^{x^ - 1) 
where 

r{p) — x^{x^ — l)Wronskain{p, x) 

does not depend from x, and our eigenvalues are zeroes of the function r{p). 

V3 — 1 

Since we can compute the Wronskian at any point 1 < a; < ^ — , using the 

V2 

series, the function r{p) is also well- determined, and we can examine where that 
function have zeroes, for examle by plotting it and looking where it have zeroes, or, 
for a more precise computation, invoking for example Newton method to find a root 
numerically. 



2.1 First few states 

The computation gives the first few zeroes at -37.169908; -81.354363; -138.473573; 
-208.859215; -292.583628; -389.671368 -500.1328 -623.97315. . . . There is an 
infinite discrete set of roots, all roots are negative. It is not difficult to get to roots 
for values of p up to about —4600. Thereafter, computation is becoming increasingly 
difficult, as the coefficients in the series grow very large, before being killed by powers 
oi x — 1 or 1/x. Also, the final product, our function r(p), becomes pretty small, with 
coefficients of orders 10^^ and results of orders 10~^, and we used tricks to go ahead. 



2.2 Orthogonality 

For the discrete set of values p„ , n = 1, 2 . . . such that r{p) — Owe have normalizable 
at 1 < X < oo wave functions {Fn{x)} . Those functions are orthogonal, which here 
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says 

/ + OC 
x^{l - x^)Fn{x)Fm{x)dx = 0, n 7^ m (8) 



3 Some pictures: r(p) for QCD4 

The graph of function r{p) up to about p — —7600 for QCD4 follows. The roots of 
this function correspond to glueball eigenstate masses. It is not difficult to get to roots 
for values of p up to about —4600. Thereafter, computation is becoming increasingly 
difficult, as for example for p around 7000 there arc terms of order lO'^"^ in the series, 
and although due to powers of of a; — 1 or l/x present the series is convergent, the 
first terms are very big, and in computing r{p) they combine into an expresson of 
order 10~^; as brut-force computation failed here, tricks were used to make it run 
for such (and higher) eigenvalues. The eigenfunctions look pretty weird here, and 
it shouldn't take much for them and eigenvalues to be corrected by about anything; 
thus we suspect our hunt for higher eigenvalues is for sport, mostly. Surprisingly, 
function r(p) keeps to be pretty nice there. ( Assistance of I.V.G. in putting up the 
pictures and fighting with latex is gratefully acknowledged.) 
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r(p) X 2 10^ 
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